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1. Introduction 

The elucidation of the universal scaling properties of equal-time and time-dependent 
correlation functions in the vicinity of a critical point was one of the most important 
achievements of statistical mechanics over the past forty years. The analogous 
systematization of the power laws and associated exponents that govern the behaviours 
of structure functions in a turbulent fluid, or in a passive-scalar advected by such a 
fluid, is a major challenge in the areas of nonequilibrium statistical mechanics, fluid 
mechanics, and nonlinear dynamics. The power-law behaviours of equal-time structure 
functions have been studied in detail over the past few decades [1] ; and, especially in the 
case of passive-scalar turbulence [2] , significant progress has been made in understanding 
the multiscaling of equal-time structure functions. The nature of multiscaling of time- 
dependent structure functions has been examined recently [3, 4, 5, 6, 7] but only for 
the case of statistically steady turbulence. We develop here the systematics of the 
multiscaling of time-dependent structure functions for the case of decaying fluid and 
passive-scalar turbulence [8]. 

To set the stage for our discussion of time-dependent structure functions in 
turbulence, it is useful to begin by recalling some well-known results from critical 
phenomena [10, 11]: At a critical point for a spin system in d dimensions, the equal-time, 
two-spin correlation function g and its spatial Fourier transform g assume the following 
power-law scaling forms: 



Here i = (|T — Tc\)/Tc., T and Tc are the temperature and the critical temperature, 
respectively, h = H/ksTc, H is the external field, ks is the Boltzmann constant, the 
spins are separated by the vector r [r = |r|], k is the wavevector, k = |k|, u, A, 
r] are critical exponents, and G and G are scaling functions. Away from the critical 
point such correlation functions decay exponentially; the associated correlation length 
diverges in the vicinity of the critical point; e.g., as ~ ii h = 0. Time-dependent 
correlation functions also assume scaling forms in the vicinity of the critical point and 
the characteristic relaxation time r diverges as suggested by the dynamic- scaling Ansatz 



which introduces the dynamic- scaling exponent z. 

The generalisation of such a dynamic-scaling Ansatz to the case of homogeneous, 
isotropic turbulence is our prime concern here. The power-law behaviours of equal- 
time structure functions, in the inertial range (to be defined later), in turbulence are 
reminiscent of the algebraic dependence on r of critical-point correlation functions. 
However, there are important differences between the two that must be appreciated 
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before we embark on a systematization of time-dependent structure functions in 
turbulence. We begin with the increments of the longitudinal component of the 
velocity 5m||(x, r,t) = [u(x + r,t) — u(x, t)] ■ (r/r) and passive-scalar 56'(x, r,i:) = 
[6'(x + r, t) — ^(x, t)], respectively; here u(x, t) and 6'(x, t) denote, respectively, the 
velocity of the fluid and the passive-scalar density at the point x and time t, and the 
subscript || the longitudinal component. The order-p, equal-time structure functions, for 
the fluid (superscript u) and passive-scalar (superscript 6) fields, are defined as follows: 

S;{r)^ ([5n||(x,r,t)f>~r^J^; 

SI{t)^ ([5^(x,r,t)n~/^ (3) 

the angular brackets indicate averages over the steady state for statistically 
steady turbulence or over statistically independent initial configurations for decaying 
turbulence; for stochastic differential equations, like the Kraichnan Model {Section 2.1), 
the angular brackets denote an average over the statistics of the noise; and the power 
laws, characterised by the equal-time exponents Cp and Cp, hold for separations r in the 
inertial range 77^ ^ r ^ L, where 77^ is the Kolmogorov dissipation scale and L the large 
length scale at which energy is injected into the system. 

Kolmogorov's phenomenological theory [1, 12, 13] of 1941 (K41) suggests simple 
scaling, with (p'^^^ = p/3, but experimental and numerical evidence favours equal-time 
multiscaling with and nonlinear, convex, monotone-increasing functions of p. For 
the simplified Kraichnan model [2, 14, 15, 16] of passive-scalar turbulence {Section 2.1) 
multiscaling of equal-time structure functions can be demonstrated analytically in 
certain limits. The analogue of the K41 theory for passive-scalar turbulence is due to 
Obukhov and Corrsin [17, 18]; if the Schmidt number 5*0 = z//k ^ 1, then their theory 
yields K41 exponents for the passive-scalar case; here u is the kinematic viscosity of the 
fluid and k is the diffusivity of the passive scalar. 

A straightforward extension of simple, K41 scaling to time-dependent structure 
functions implies that the dynamic exponents z^'^^ = 2/3 for all p. This naive extension 
fails for two reasons: (a) it does not distinguish between the temporal behaviours of 
structure functions of Eulerian, Lagrangian, and quasi-Lagrangian {Section 2) velocities 
or passive-scalar densities; and (b) it does not account for the multiscaling of structure 
functions. These difficulties have been overcome to a large extent for statistically 
steady turbulence [3, 4, 5, 6, 7, 19] as we summarise below. There is consensus now 
that Eulerian structure functions display simple scaling with only one dynamic-scaling 
exponent = 1 because of the sweeping effect: the mean flow, or the flow caused by 
the largest eddy, advects small eddies, so spatial separations r in (3) are related linearly 
to temporal separations r via the mean-flow velocity [19]. By contrast, it is expected 
that Lagrangian [19] or quasi-Lagrangian [3, 4, 5, 6, 7, 20] time-dependent structure 
functions should show nontrivial dynamic multiscaling. The task of extracting well- 
averaged time- dependent Lagrangian or quasi-Lagrangian structure functions from a 
direct numerical simulation (DNS) of the Navier-Stokes equation is a daunting one [21]: 
a dynamic exponent has been extracted from a full Lagrangian study [19] only for order 
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p = 2. Thus the elucidation of dynamic multiscahng has rehed on predictions based on 
generahsations of the multifractal formahsm [3, 4, 5, 6, 7] and on numerical studies of 
shell models [4, 5, 6, 7]. These studies show that, if dynamic multiscaling exists, time- 
dependent structure functions must be characterised by an infinity of time scales and 
associated dynamic multiscaling exponents [6]. Furthermore, the dynamic exponents 
depend on how we extract time scales from time-dependent structure functions; e.g., for 
fluid turbulence, time scales obtained from integrals (superscript / and subscript 1) and 
second derivatives (superscript D and subscript 2) of order-p time-dependent structure 
functions yield the different dynamic exponents z^^ and ■ Finally, the different 
dynamic multiscaling exponents are related by different classes of linear bridge relations 
to the equal-time multiscaling exponents. For a careful discussion of these issues we must 
of course define time-dependent structure functions. The details necessary for this paper 
are given in Section 2 and Section 4. 

The dynamic multiscaling of time-dependent structure functions described briefly 
above applies to statistically steady turbulence. Does it have an analogue in the case 
of decaying turbulence, since time-dependent structure functions must, in this case, 
depend on the origin of time to at which we start our measurements? This question 
has not been addressed hitherto. We show here how to answer it in decaying fluid and 
passive-scalar turbulence [8]. In particular, we propose suitable normalisations of time- 
dependent structure functions that eliminate their dependence on to! we demonstrate 
this analytically for the Kraichnan version of the passive-scalar problem and its shell- 
model analogue and numerically for the GOY shell model [1, 22, 23] for fluids and a 
shell-model version of the advection-diffusion equation. In these models we then analyse 
the normalised time-dependent structure functions for the case of decaying turbulence 
like their statistically steady counterparts [6, 7]. This requires a generalisation of the 
multifractal formalism [1] that finally yields the same bridge relations between dynamic 
and equal-time multiscaling exponents as for statistically steady turbulence [6, 7]. For 
the Kraichnan version of the passive-scalar problem we show analytically that simple 
dynamic scaling is obtained. This is because [Section 3) the advecting velocity is 
random and white in time. In addition, we find numerically for shell models of fluid and 
passive-scalar turbulence that dynamic-multiscaling exponents have the same values for 
both statistically steady and decaying turbulence; so, in this sense, we have universality 
of the multiscaling of time-dependent structure functions in turbulence. The equal-time 
analogue of this universality has been discussed in Ref . [24] . 

The remaining part of this paper is organized as follows. In Section 2 we introduce 
the models we use and give the details of our numerical simulations. Section 3 presents 
our analytical studies of decaying turbulence in the Kraichnan model and its shell-model 
analogue. Section 4 shows how to generalise the multifractal formalism to allow for time- 
dependent structure functions in decaying turbulence and how to obtain bridge relations 
between dynamic and equal-time multiscaling exponents in this case. In Section 5 we 
present the results of our numerical studies of dynamic multiscaling in the GOY shell 
model for fluid turbulence and for shell models of a passive-scalar field advected by a 
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turbulent velocity field. Section 6 is devoted to a discussion of our results in the context 
of earlier studies; we also suggest possible experimental tests of our predictions. 

2. Models and Numerical Simulations 

We have used several models to study time-dependent structure functions in fluid and 
passive-scalar turbulence. These range from the Navier-Stokes and advection-diffusion 
equations to simple shell models; the latter are well-suited for our extensive numerical 
studies. It is useful to begin with a systematic description of these models. 

Fluid flows are governed by the Navier-Stokes (NS) equation (4) for the velocity 
field u(x, t) at point x and time t, augmented by the incompressibility constraint (5), 
since we restrict ourselves to low Mach numbers: 

(9tU + u.Vu = -VP + i/oV^u + f; (4) 

V.u = 0. (5) 

Here vq is the kinematic viscosity, P the pressure, the density p is taken to be 1, and 
f the external force, which is absent when we consider decaying turbulence. If i and 
V are, respectively, characteristic length and velocity scales of the flow, the Reynolds 
number Re = ^ provides a dimensionless measure of the strength of the nonlinear term 
in (4) relative to the viscous term; for the case of decaying turbulence it is convenient to 
use the Reynolds number for the initial state, i.e.. Re with v the root-mean-square 
(rms) velocity of the initial condition and i the system size [the linear size of the 
simulation box in a direct numerical simulation (DNS)]. Given the incompressibility 
condition (5), the pressure can be eliminated from (4) and related to the velocity by a 
Poisson equation. The equation for the velocity alone is most easily written in terms of 
the spatial Fourier transform u(k, t) of u(x, t); and it can be shown easily that u(k, t) 
is affected directly by all other Fourier modes. This is the mathematical representation 
of the sweeping effect in which the largest eddies (i.e., modes with small /c =| k |) 
directly advect the smallest eddies (i.e., large-A; modes); such direct sweeping lies at the 
heart of the Taylor hypothesis [25] and leads eventually to trivial dynamic scaling for 
time-dependent structure functions of Eulerian fields with a dynamic exponent = 1 
for fluid turbulence. 

As we have mentioned above, nontrivial dynamic multiscaling is expected if we use 
Lagrangian or quasi-Lagrangian velocities. The Lagrangian formulation is well known 
[26]; the quasi-Lagrangian [3, 20] one uses the following transformation for any Eulerian 
field V(x, t): 

^(x,t) =V^[x + R(t;ro,0),t], (6) 

where is the quasi-Lagrangian field and R(t; Tq, 0) is the position at time t of a 
Lagrangian particle that was at point Tq at time t = 0. 
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The advection-diffusion (AD) equation for the Eulerian passive-scalar field 6{'x^t) 

is 

— + u.ve = KV'e + (7) 

where k is the passive-scalar diffusivity and, if we consider decaying passive-scalar 
turbulence, the external force fg is set to zero. The advecting velocity field u should be 
obtained, in principle, by solving equations (4) and (5). By using equation (6) we get 
the quasi-Lagrangian version of the advection-diffusion equation (7): 

+ [u(x, t) - u(x, 0)] . V^(x, t) = kS/H{^, t) + /,(x, t). (8) 

Direct numerical simulations of equations (4) and (5) or equation (7), though 
feasible, have not yet provided data that are averaged well enough to yield reliable time- 
dependent structure functions of quasi-Lagrangian velocity [21] or passive-scalar fields. 
Time-dependent Lagrangian structure functions have been obtained [19] only for order 
p = 2. Thus a first-principles DNS study of dynamic multiscaling in fiuid or passive- 
scalar turbulence is not possible at the moment. However, significant progress has been 
made in statistically steady turbulence by studying dynamic multiscaling in simplified 
models like the Kraichnan model for passive-scalar turbulence and shell models for 
fluid and passive-scalar turbulence. We discuss these models below since our studies of 
decaying turbulence will be based on them. 

2.1. The Kraichnan Model (Model A) 

The Kraichnan model for passive-scalar turbulence [2, 14, 15, 16] begins with the 
advection-diffusion equation (7) but replaces the Navier-Stokes velocity field by one in 
which each component Mj(x, t) of the velocity is a zero-mean, delta-correlated, Gaussian 
random variable with the covariance 

(Mi(x,t)Mj(x + r,t')) =2Dij{Y)5{t-t'). (9) 



The Fourier transform of Dij{r) has the form 

2 



/ 1 \ -(.d+t,)/^ 



A -Mi 

9 



(10) 



where q is the wave vector, L the characteristic large length scale, rja the dissipation 
scale, and ^ a tunable parameter. In the limits L — > oo and 77^ ^ 0, of relevance to 
turbulence, we have, in real space. 



D,,{r)=D%,,-U,{r), (11) 



2 

with 

di^ = Dir^ 



{d-l+ 05., - , (12) 

where ~ CiL^; Ci and Di are dimensional constants. We refer to (7-12) as Model 
A to distinguish it from other models that we use. For < ^ < 2, this model shows 
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multiscaling of order-p, equal-time passive-scalar structure functions, as can be shown 
analytically, in certain limits [2]. However, for the case of statistically steady turbulence, 
this model exhibits simple dynamic scaling [6]. 

2.2. Shell Models 

We will also use some shell models for fluid and passive-scalar turbulence. These 
models are highly simplified representations of the Navier-Stokes or the advection- 
diffusion equations (4-7) and are, therefore, far more tractable numerically than (4-7). 
Nevertheless, shell models retain enough properties of their parent equations to make 
them useful testing grounds for the multiscaling of structure functions in turbulence. 
Shell models are defined on a logarithmically discretised Fourier space in which complex 
scalar variables (e.g., the velocity n„ or passive scalar 9^) are associated with the shells 
n and scalar wave vectors = /cqA"; typically A = 2, fco = 1/16; and the boundary 
conditions are that the shell variables vanish if n < 1 or n > N (we use = 22). 
These models consist of coupled, nonlinear, ordinary differential equations (ODEs) that 
specify the temporal evolution of the shell variables Un and 6n- Shell-model ODEs are 
similar to the Fourier-space versions of their parent partial differential equations: (a) 
their dissipative terms are linear in one of the shell variables and quadratic in /c„; (b) 
their analogues of advection terms are linear in kn and bilinear in the shell variables; 
e.g., for fluid turbulence a representative term is of the form i/c„M„M„/, with n ^ n'; 
and (c) they conserve the shell- model analogues of the energy, helicity, etc., in the 
absence of dissipation and forcing. However, variables in a given shell are influenced 
directly only by their nearest- and next-nearest-neighbour shell variables; by contrast, 
Fourier transformations of the NS and the AD equations couple every Fourier mode to 
every other Fourier mode, leading to the sweeping effect mentioned above. Thus direct 
sweeping is absent in shell models, so they are often thought of as an approximate, 
quasi-Lagrangian representation of their parent equations. 

For studies of decaying turbulence one can envisage several initial conditions. We 
have used initial conditions of two types: (a) in the first (Type-I) we drive the system 
to a statistically steady turbulent state by forcing the first shell {n = 1); we then 
turn off the force and allow the turbulent state and the associated energy spectrum to 
decay freely; our measurements are made in this decaying state; and (b) for the case of 
fluid turbulence we use a second initial condition (Type-II) in which all the energy is 
concentrated in the first few shells with small kn, i.e., large length scales; we then allow 
the system to evolve without any force; the energy cascades to large values of /c„ till 
the energy spectrum becomes similar to that in forced turbulence; this spectrum then 
decays slowly in time and the measurements we report are made during this stage of 
evolution. 
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2.3. Model B 

We use the shell-model analogue of the Kraichnan model introduced in [27] in which 
the equation for the passive-scalar variable On is 



d 
It 



6„ 



+ K 



'n+2'^n+l 



+ fn 



(13) 



where the asterisks denote complex conjugation, a„ = /Cn/2, 6„ = — A;„_i/2, and 
Cn = kn+i/2; fn is an additive force that is used to drive the system to a steady state; the 
boundary conditions are m_i = Uq = = = 0',un+i = = ^tv+i = = 0. 

The advecting velocity variables are taken to be zero-mean, white-in-time, Gaussian 
random complex variables with covariance 

{Un{t)ul{t')) = C^K^SmJit - t'), (14) 

where C2 is a dimensional constant. We refer to equations (13-14) as Model B. 

In our numerical simulations of this model we first obtain a statistically steady 
turbulent state by forcing the first shell with a random, Gaussian, white-in-time force. 
The force is then switched off and measurements are made as the turbulence decays. 
We use a weak, order-one, Euler scheme to integrate the resulting Ito form [27] of (13) 
with an integration time step 6t = 2~^^, diffusivity k = 2~^^, and ^ = 0.6. 

For such a passive-scalar shell model the order-p, equal-time, structure function 
and its exponent are defined via 

s'^{K)^{%it)e:{t)r/')r^k-'^'; (15) 

it is natural, therefore, to define the time-dependent version of Sp{kn) as follows: 

F^ik^,to,t) ^ {[9nito)0:ito + t)r/'). (16) 

The power-law dependence on the right-hand-side of (15) is obtained for kn in the inertial 
range. In our numerical calculations we use extended self-similarity (ESS) to extract 
the exponent ratios Cp/Ci- 



2.4. Model C 

The most commonly used shell-model analogue of the NS equation is the GOY 
model [1, 22, 23]: 



The coefficients an 



Ur 



(17) 



kni bn 



-5kn- 



5)kn-2 are chosen in a manner 



that conserves the shell-model analogues of energy and helicity in the inviscid, unforced 
limit; an external force fn drives the system to a steady state. We use the standard 
value 5 = 1/2; the boundary conditions are u_i = Uq = 0]Un+i = = 0. We will 

refer to this as Model C. 
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We use two different kinds of initial conditions in our study of decaying fluid 
turbulence in this model. For Type-I initial conditions we first drive the system to a 
statistically steady turbulent state with an external force /„ = x 5 x 10~^5„,i. The 

force is then switched off and the shell velocities at this instant are taken as the initial 
condition. The turbulence then decays. Our structure-function measurements are made 
during this period of decay. To obtain the second type of initial condition, the energy is 
initially concentrated in the first few shells by choosing the following initial (superscript 
0) velocities: m° = kl/'^e''^", for n = 1,2, and m° = kl/'^e'^"^ e'-'^" , ior 3 <n < N, with 
"dn a random phase angle distributed uniformly between and 27r. This energy then 
cascades down the inertial-range scales without significant dissipation until it reaches 
dissipation-range scales at cascade completion. The energy dissipation-rate per unit 
mass shows a peak, as a function of time, roughly at cascade completion [9, 28], and 
the energy spectrum E{k) and the structure function (18) show well-developed inertial 
ranges. These decay very slowly in time, so at each instant exponents can be determined 
from plots of E{k) and the structure functions. Therefore, for initial conditions of this 
type, we wait for cascade completion before making measurements of structure functions. 

We employ the slaved, Adams-Bashforth scheme [29, 30] to integrate the GOY- 
model equations with a time step St = 10~^. In our numerical simulations the viscosity 
u = 10~^ and the total number of shells = 22; this provides us with a large inertial 
range from which exponents can be obtained reliably. 

For this model, the order-p, equal-time structure function and its exponent are 
defined as follows: 

S;ik^) ^ (K(t)<(t)]*'/2) ~ kn^^; (18) 

the associated time- dependent structure function is 

F;{kn,to,t) = (K(to)<(to + t)]^/'>, (19) 

where the power-law dependence on the right-hand-side of (18) holds for kn in the inertial 
range. (For statistically steady turbulence, the time-dependent structure function has 
no dependence on to, so, without loss of generality, to can be taken to be 0.) A direct 
determination of from (18) is not very accurate because of an underlying 3-cycle in 
the static version of the GOY shell model [31]. The effects of this 3-cycle can be filtered 
out to a large extent by using the modified structure function 

S^(n) = {\^[u^+2Un+iUn ' (1/4)m„„iM„m„+i] ^/^^ ~ k~^^ ; (20) 

we use Sp(n) in our numerical calculation of (p. We measure time in terms of 
the initial large eddy-turnover time ti = 1/ (urmski); the root-mean-square velocity 

2.5. Model D 

A turbulent velocity field does not have the simple statistical properties assumed in 
Models A and B. To overcome this we study the shell model of Ref. [32], hereafter 
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referred to as Model D, in which the advecting velocity Uyi IS Cl solution of the GOY 
shell model(17). The passive-scalar shell variables 9^ obey 



d 2 




Jt + ""^^ 


On = i 



an{9, 



-1 "n—1 n {On— \Un-2 + On-2Un-l) 



+ Cn{9n+2Un+l + 9n+lUn+2) + /) 



(21) 



where a„ = /c„, hn = -kn-i/2, and c„ = -kn+i/2; /„ = (1 + «) x 5 x 10 ^6n,i is an 
additive force that drives the system to a steady state; the boundary conditions are 
U-i = Uo = 9_i = 9q = 0; ujv+i = njv+2 = ^'at+i = 9N_^.2 = 0. 

For this model we start with Type-I initial conditions, i.e., we force both the coupled 
equations (17) and (21) till a statistically steady turbulent state is obtained and then 
switch off the force. The shell variables at this instant of time are taken as the initial 
condition; and then the turbulence is allowed to decay. 

We employ a second-order Adams-Bashforth scheme to integrate Model D with 
a time step 6t = 10~^ and set the diffusivity k = 5 x 10~^ so the Schmidt number 
u/k, = 1/5; and = 22 as in Model C. The definitions of structure functions for Model 
D are the same as those for Model B, i.e., equations (15) and (16). We limit the effects 
of the 3-cycle (mentioned above for Model C) in our numerical evaluations of (p by using 
the modified structure function 



S^(n) = {\^[9n+29n+l9n " {I / A)9n-l9n9n+l. 



(22) 



3. Analytical Results for Models A and B 

In this Section we present our analytical results for time-dependent, passive-scalar 
structure functions for the Kraichnan model (Model A) and its shell-model analogue 
Model B. For Model A we obtain results for both Eulerian and quasi-Lagrangian 
structure functions. We find, in particular, that time-dependent structure functions 
for these models can be factorised into a part that depends on the time origin to and 
a part that depends on t but is independent of to- This important result motivates a 
similar factorisation hypothesis that we propose, and verify numerically, for Models C 
and D in subsequent Sections. 



3.1. Model A 

Consider first the Eulerian version of the AD equation (7). We assume that a turbulent 
statistical steady state has been established because of an external force fe. We turn 
off this force at time 0. A spatial Fourier transform of (7) now yields 

^^^^ = I j %Sj(q, t)^(k - q, t) A - Kkjkj9{k, t), (23) 

where the tildes denote spatial Fourier transforms, we sum over repeated indices, and 
the statistics of Uj{q,t) are specified by equations (9) and (10). The Fourier-space, 
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second-order correlation function jF|(k, to,t) = (6'(— k, to)^(k, to + t)) {to > is any- 
time origin) satisfies the equation 

^-^l(k,to,t) ^ ^ ^ ^ de{k,to + t) ^ 

which can be combined with (23) to get 

^•^'^^^^"'^^ = tk, j {e{-k, to)uj{q, to + t)e{k - q, to + t))d'q 

- Kkjkj{e{-k,to)9{k,to + t)). (25) 

We average over the statistics of the advecting velocity field (9) by using Novikov's 
theorem [33]: e.g., the first term in (25) reduces to 



(^(-k, to)n,(q, to + t)^(k - q, to + t)) 



dt' 



(M,(q, to + t)Mi(-q, t')) (^~(-k, to)—- — — 7^(k - q, to + t')) 

5ui(-q,to + t') 



(26) 



Finally equations (9) and (23)-(24) yield 



dTlik, to, t) ^ _2fc,fc,^2'(k, to, t) r AA,(g). (27) 



Since 2 Dijd'^q = D^6ij ~ CiL^, where Ci is a dimensional constant, in the limits 
ft — > 0, r/rf — > 0, and L ^ oo of relevance to turbulence, we get 

dJ'^ir^to^t) .d^fl^ryto^ 

A spatial Fourier transform allows us to integrate this equation to obtain 

:F^{k, to, t) ~ ^l{k, to) exp [-CiL^kH] . (29) 

If we set t = 0, we see that jF|(k, to,0) is just the passive-scalar, equal-time structure 
function; ip2{k,to) is, therefore, proportional to the Fourier transform of the passive- 
scalar equal-time structure function S2{r). Thus, for a fixed but large value of L, we 
get, in the Eulerian framework, simple dynamic scaling with an exponent zf = 2. We 
note that we get an exponent which is not equal to unity {Section 1) in this case because 
of the white-in-time nature of the advecting field. 

The analogue of equation (28) for the second-order, quasi-Lagrangian structure 
function follows from [7] the quasi-Lagrangian form of the advection-diffusion equation 
(8). 

g^|(r,to,t) _ ,^0^ _ , a^|(r,to,t) 
dt dr.dr, 

d4{r,to,t) 

By substituting for dij from (12) we obtain, for the isotropic case, 
dH{r,to,t) ^ d'j^l{r,to,t) 

di • ^^^^ 
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A spatial Fourier transformation allows us to integrate this equation to obtain, in the 
hmits K — * 0, ?7 — > 0, and L oo, 

^2^(k, to, ~ vi{.k, to) exp [-Dip-H] . (32) 

ip2{k,tQ) is now proportional to the Fourier transform of the equal-time, quasi- 
Lagrangian, passive-scalar structure function (which is, of course, the same as the 
equal-time, Eulerian, passive-scalar structure function). Equation (32) shows that, in 
the quasi-Lagrangian framework, jF|(k, to,^) factorises into a part that depends on to 
and another which depends only on t. From the second factor we get simple dynamic 
scaling with an exponent Z2 = 2—^, which is different from the Eulerian exponent 2;f = 2 
for Model A. Such a factorisation should also follow for higher-order, time- dependent, 
structure functions as we show in the next subsection for Model B. 



3.2. Models 



We can use the methods of the previous subsection to obtain analytical expressions 
for time-dependent structure functions for Model B. Consider first the second-order 
structure function 



whence 



F^in,to,t) = {9n{to)e:ito + t)), 



- = {(^n[to) tt: /, 



(33) 



(34) 



dt dt 
the angular brackets denote an average over the statistics of u„(t) that are specified by 
equation (14). By using the complex conjugate of (13) in (34) and Novikov's theorem 
we get terms of the form 



Uto)0:+,{to + t)<_i(to + t)) = «_i(to + t)Un-l{to + t)) X 

/ 5 

-0nito)6n+l{to + t) 

Finally by using (14) we obtain 

gF|(n,to,t) ^ 1 
dt 4 

where C2 is a dimensional constant. Integration now yields 

1 



\5Un~l{tQ + 1) 

C2kl-^A{OF^{n,to,t), 



F^in.tQ.t) = 02(n,to)exp 



4 



C2kl-^A{0t 



(35) 



(36) 



(37) 



with A(0 = (2(2?-2) + 2-(2C-2)) + (2? + 2-«) + (2(«-2) ^ 2-«-2))^ Similarly, we obtain 
the following exact expression for the fourth-order structure function: 

1 



F^{n, to, t) = (j)l{n, to) exp 



'-C2kl-^A{0t 



(38) 



(j)2{n,to) and ^4(^,^0) are, respectively, the second- and fourth-order, equal-time, quasi- 
Lagrangian, passive-scalar structure functions. Thus 22 = -^4 = 2 — ^, as for the quasi- 
Lagrangian structure functions of Model A. (Recall that we expect quasi-Lagrangian 
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behaviour for shell models since they do not have a direct sweeping effect.) The equality 
of Z2 and indicates that we have simple dynamic scaling in Model B. We expect 
all the quasi-Lagrangian exponents 2;^ to be 2 — ^ for this model, but the analytical 
demonstration of this result becomes more and more complex with increasing p. 

Given a factorisation of the form shown in (37), it is possible to normalise the time- 
dependent, passive-scalar structure function Fp{n,to,t) by its value at t = and thus 
make it independent of to- We cannot prove that such a factorisation exists in Models 
C and D, but we present compelling numerical evidence for it in Section 5. 



4. Multifractal Formalism for Models C and D 



Equal-time Eulerian and quasi-Lagrangian structure functions are the same for 
homogeneous and isotropic turbulence [34]. Since quasi-Lagrangian structure functions 
are required for our study of dynamic multiscaling, we present the multifractal formalism 
in terms of quasi-Lagrangian variables. Multiscaling in statistically steady fluid 
turbulence can be rationalised by using the multifractal formalism, which assumes that 
a turbulent flow has a continuous set of scaling exponents h in the set I = {hmini hmax), 
instead of a single exponent (e.g., h = 1/3 yields simple K41 scaling) [1]. The scaling 
exponents h characterise the behaviour of velocity differences ^^^(x): For each /i G X 
there exists a set S/^ C of fractal dimension T>^{h) and 5ur{^)/uL ~ {r / L)^ for 
separations r in the inertial range if x e Sh and with ul the velocity at the forcing scale 
L. Given the measure dn{h) for the weights of the different values of h, the order-p, 
equal-time velocity structure function is 

ui ul J J L 

where the ph term comes from p factors of (r/L)^ and the additional factor of 
(r/L)^"-^"'^'^) is the probability of being within a distance ~ r of the set H^, of dimension 
V^{h), which is embedded in three dimensions. V^[h), hmin, and hmax are assumed to 
be universal. In the limit r/L 0, of relevance to fully developed turbulence, we get 
the equal-time scaling exponent = mfh[ph + 3 — V^{h)] by the method of steepest 
descents. 

We now define the order-p, time-dependent, quasi-Lagrangian velocity structure 
function 

^;(r,{ti,...,tp}) = ([5M||(x,r,ti)...5M||(x,r,tp)]). (40) 

For simplicity, we consider ti = t and t2 = ■ ■ ■ = tp = 0, denote the structure function by 
J-'p{r,t), and suppress the subscript ||, i.e., we use ^^^(x) = 6u\\{'x,r,ti) . The natural 
extension of the multifractal formalism to the case of time-dependent structure functions 
in statistically steady fluid turbulence follows from the Ansatz [6] 

««/...(.)— aW^), 
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where the scahng function ^P''* (^t^^ is assumed to have a characteristic decay time 

Tp^h ~ r/5M^(x) ~ r^~^ and ^^'^(0) = 1. 

For statistically steady passive-scalar turbulence the application of the multifractal 
formalism is more complicated than it is for fluid turbulence. This is because we must 
now deal with a joint multifractal distribution of both the velocity and passive-scalar 
variables. Therefore, the order- equal-time structure function for passive-scalar 
turbulence 

STAt) ^ (5^f(x)5<'(x)) (42) 
has the multifractal representation 

MK9)C-,T'''''-^'''''-'\ (43) 



where u and 9 are assumed to possess a range of universal scaling exponents /i G X = 
{hmin, hmax) and g eX' = {grain, gmax), respectively. For each pair of h and g in these 
ranges, there exists a set "Sh^g C M.^ of fractal dimension V^'^{h,g). The increments 
in the velocity 5ur(x) and the passive-scalar field 66r{'x.) scale as 6{ir{'x.)/uL ~ {r / L)^ 
and 66r{'x)/6L ~ {r/L)^ for separations r in the inertial range if x G S/j^^. Hl and 6l 
are, respectively, the velocity and the passive-scalar variables at the forcing scale L. 
For simplicity we will only consider passive-scalar structure functions with = in 
equation (42), i.e., 

S'pir) = (5^J?(x)), (44) 

which have the multifractal representation 

V / dfxih,g)(-) ■ (45) 



§1 Jxj, \LJ 

as before, the equal-time exponents (p can be related to V^'^{h,g) by the method of 
steepest descents. The corresponding order-p, time-dependent passive-scalar structure 
function is 

J^'p{r,{tu...,tp}) = {m^,r,t,)...69{^,r,tp)]). (46) 

As in (40), we consider ti = t and t2 = ■ ■ ■ = tp = for simplicity. Given 
the nature of multiscaling in passive scalars advected by a turbulent velocity field, 
it is possible to understand passive-scalar turbulence within the framework of the 
multifractal formalism. However, the analogous expression for time-dependent structure 
functions in passive-scalar turbulence has to take into account the multifractal nature 
of the advecting velocity field [7]. We generalise the multifractal representation of time- 
dependent structure functions in the following way: 

oc / di,{K 9) (j) G'''^' ^ ; 47 

9^r J 1,1' ^LJ \rphgj 
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and we assume that (a) the function QP'^'9(^—L- 

'''phg 

(b) Q^'^'^{0) = 1, and (c) that the dominant contribution to r, 



has a characteristic decay time Tphg, 
phg, in the hmits k ^ 0, 



r/d — > and L ^ oo, has the following scaling form 



Tphg ~ r/5n,(x) ~ r^-^ (48) 

Given Tp{r,t) and J-'^{r,t), we define the order-p, degree-M, integral- (/) and 
derivative- (Z)) times as follows [6, 7] (if the integrals and derivatives in equations (49) 
and (50) exist): 

(l/Af) 



and 



rp%%r) 



St{r) Jo 

1 d^'T^{r,t) 



J^i{r,t)t^^-''>dt 



dt^^ 



t=o 



(-1/A/) 



(49) 



(50) 



exponents for fluid turbulence zi'^ and zfii" via 



Tp'^ir) 



I,u 



respectively, with either u or 6. For statistically steady turbulence we can then use 
the dynamic-multiscaling Ansatz to define the integral- and derivative-time multiscaling 

p,M ^^'^ ^p,M 

(51) 
and 

^,5/(0 - r<i'. (52) 

respectively, for r in the inertial range. 

By substituting the multifractal form (41) in (51), evaluating the time integral first, 
and then performing the integration over the multifractal measure by the saddle-point 
method, we obtain the integral bridge relations 



I,u 



1 + [Cm - C] /M, 



(53) 



which were first obtained in [3]. Similarly we get the derivative bridge relations 



-p,M - 1 + [C; - C+m] /M, (54) 

which were first obtained in forced Burgers turbulence [35, 36] for the special cases (a) 
p = 2, M = 1 and (b) p = 2, M = 2, respectively, and for statistically steady fluid 
turbulence in [6]. We note that within the K41 phenomenology (Cp'^'^^ = p/3) the 
bridge relations yield the same dynamic exponent z^^^ = 2/3 for both integral- and 
derivative-time scales. 

For passive-scalars advected by a turbulent velocity field, the corresponding 
dynamic-multiscaling exponents are defined via 



(55) 



and 



(56) 
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To obtain bridge relations for dynamic-multiscaling exponents in statistically steady 
passive-scalar turbulence we define the degree-M, order-p integral-time for this case: 



dfi{h,g) 



1,1' 



L 



M 

p,g,h 



-',P,h,g 



t 



Tpgh 



d 



t 



Tpgh 



1 



(57) 



where Z = ?> + pg — V^'^{h,g) and the argument of the scaling function is suppressed 
for notational convenience. By using the scaling form of Tpgh (48) and assuming, as in 
Ref.[32], that 

{6eP6u-'') ^ {6eP){6u-''), (58) 

we get 



l/M 



(59) 



/ 6 

Thus the degree-M, order-p, integral-time, dynamic-multiscaling exponent Zp\ = 
1 + Cm/M. 

Similarly we obtain the degree-M, order-p derivative time scale for passive-scalars 
by substituting for JF^ from (47) in (50) and using (48) and (58): 

-1 -l/M 



■^p,M V ) 




1 



1 



l~Qj/M 



dfi{h,g) 
I d^i{h,g){^ 



d 



M 



dt^^ 

dt^^ 

-1 -l/M 



gp,h,g 
gp,h<g 



1 



t=0/ T, 



t=0 



pgh^ 



-1 -l/M 



-l/M 



(60) 



which, along with (56), yields the bridge relation zf*'' 



^pM ~ ^ ~ Cm/M. To summarise, 
the bridge-relations for degree-M, order-p, derivative- and integral-time dynamic- 
multiscaling exponents are, respectively: 

i,e 1 , <5-M 



D.6 



/■u 

M' 



-p,m - m ' ^P''' ~ ' M' ^^^^ 
Note that, in contrast to the bridge relations (53) and (54), there is no p-dependence 
on the right-hand sides of these relations. However the M-dependence is the signature 
of nontrivial dynamic multiscaling here. Furthermore, the integral scale bridge relation 
is meaningful only for those values of M for which (^j^^j is well defined [7]. We now 
extend our discussion to the the case of decaying turbulence and define the order-p. 
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time-dependent quasi-Lagrangian velocity structure function (for fluid turbulence) and 
passive-scalar structure function (for passive-scalar turbulence) as 

^;(r,{ti,...,tp}) = ([5n||(x,r,ti)...5n||(x,r,g]) (62) 

and 

^^(r, {ti, . . . , t^}) = ([(5^(x, r, t,)...69{^, r, tp)]), (63) 

respectively. For simplicity, we consider the case ti = to + t and t2 = . . . = tp = t^. 
The derivations of the bridge relations we have given above go through if we assume 
the following multifractal forms for time-dependent velocity and passive-scalar structure 
functions, respectively, in decaying turbulence: 

oc A%r, <„) / Mn) (9'^"-"'" a- (^ I ; (64) 

and 



i.e., we assume that the multifractal form factorises into a part A"(r, to) [or A^(r, to) in 
the case of passive scalars], which depends on the origin of time to, and an integral that 
is independent of to. This assumption is motivated by the factorisation we have seen in 
Section 3 above in (32), (37), and (38) for time-dependent passive-scalar structure 
functions in Models A and B. [The analogous expressions for equal-time structure 
functions in decaying turbulence are obtained by setting t = in (64) and (65).] Given 
(64) and (65), the integral- and derivative-time scales for decaying turbulence become 
independent of to and assume forms identical to (49) and (50); consequently the bridge 
relations (53), (54) and (61) remain unaltered in decaying turbulence. 

In the next Section we give compelling numerical evidence in support of our 
assumptions (64) and (65). 



5. Numerical Results 



We now present our numerical results for djTiamic multiscaling in decaying, 
homogeneous, isotropic turbulence in Models B,C, and D in Section 5.1, Section 5.2, 
and Section 5.3, respectively. We show en passant that, for Models B and D, equal- 
time multiscaling exponents are universal in the sense that the exponents obtained from 
decaying-turbulence runs are equal (within error bars) to the exponents for statistically 
steady turbulence in these models [6, 27, 32]. The underlying reason for this universality 
for passive-scalar turbulence has been discussed earlier [37, 38], but we believe that ours 
is the first numerical demonstration of such universality in these models. Since we use 
shell models in this Section, we replace u and 6' by m and 9, respectively. 
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5.1. Model B 

Let us begin with the equal-time exponents Cp for the case of decaying turbulence in 
Model B. As in Ref. [27], we use the extended-self-similarity (ESS) procedure to obtain 
exponent ratios from log-log plots [Fig. (la)] of the equal-time structure functions 
Sp{kn) versus S'2(/cn)- In particular, the slopes of the linear regions of the plots in 
Fig. (la), with 4 < n < 12, yield the inertial-range exponent ratios Cp/Cl- From 50 
such statistically independent runs (each run, in turn, is averaged over 5000 independent 
initial conditions) we calculate the means of the equal-time multiscaling exponent ratios 
Cp/C2j these are plotted in Fig. (lb) as functions of the order p for 1 < p < 8. The 
standard deviations of these exponent ratios (calculated from our 50 runs) provide 
us with error bars; these are smaller than the symbol sizes used in Fig. (lb). Our 
exponents ratios are in agreement (within the error bars) with earlier results for equal- 
time exponents for statistically steady passive-scalar turbulence in Model B [7, 27]. 
We have also determined the exponent (2 directly from log- log plots of T,2{kn) versus 
kn [equation (22)]. We find C2 = 1-403 ± 0.003; this agrees well with the analytical 
prediction [27, 39] = 2 — since we use C, = 0.6 in our simulations. 




Figure 1. (a) Log-log plots of the order-p structure functions Sp{kn) versus the 
second-order structure function S2(kn) from our numerical simulations of Model B for 
p = 1 (uppermost curve) to p = 8 (lowermost curve) . (b) Plot of the exponent ratios 
Cp/C2i obtained from the slopes of the linear, inertial region in (a), versus p. Our data 
points (o) are connected by a line to guide the eye; the error bars are smaller than the 
sizes of our symbols; the dashed straight line corresponds to the analogue of the K41 
scaling prediction for this case. 

We begin the discussion of our numerical results for time-dependent, passive-scalar 
structure functions in decaying turbulence by comparing our analytical result (37) with 
data from our simulations of Model B. In Fig. (2a) we show the agreement between the 
two. The error bars have been calculated as in the the case of equal-time multiscaling 
exponents: The data we present are the means of the values obtained from 50 different 
statistically independent runs; and the standard deviations of these values yield the error 
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bars. Here and henceforth time is scaled by the large-eddy-turnover time t^. Given the 
exponential decays of time- dependent passive-scalar structure functions in Models A 
and B [equations (32) and (37)], we can extract a unique time scale from plots like 
the one in Fig. (2b). Thus we do not expect time-dependent structure functions to 
exhibit dynamic multiscaling in Models A and B. In particular, equation (37) implies 
Tiikn) = [l/4/c^"'^y4(,^)]"^ whence Z2 = 2 — ^. For our choice of ^ = 0.6 we should, 
therefore, obtain Z2 = 1.4, and indeed our numerical simulations yield Z2 = 1.398±0.003 
[Fig. (2b)]. To demonstrate numerically that higher-order time-dependent structure 
functions also have a dynamic scaling exponent equal to Z2, as shown in Section 3.2 
for the fourth-order time- dependent structure function, we show in Figs. (3a) and (3b) 
our numerical results for the third- and fourth-order time-dependent structure functions 
versus t/t^ for 6 < n < 13 from which we extract T^{kn) and T^{kn), respectively. The 
insets show log- log plots of these times versus kn, the slopes of these plots yield the 
dynamic-scaling exponents z^ = 1.398 ±0.003 and z^ = 1.402 ±0.005 in agreement with 
our expectation Zp = 2 — ^, for all p, since ^ = 0.6 in our calculations. 




Figure 2. (a) Plots comparing tlic normalised, second-order, time-dependent structm'e 
function versus the dimensionless time 1/1^ from our numerical simulations (o) of 
Model B with our analytical expression (37) (*). For clarity we show data for four 
shells, from n = 6 (uppermost) to n = 9 (lowermost). The error bars on our numerical 
simulations are shown as vertical lines on the data points. 

(b) A log-log plot of T2 (fc„) versus The slope of this plot gives the dynamic scaling 
exponent Z2 = 1.398 ± 0.003. Our result from numerical simulations agree well with 
the analytical result Z2 — 2 — since we have chosen ^ — 0.6. 



5.2. Model C 

The equahty of equal-time exponents for decaying and statistically steady fluid 
turbulence was demonstrated numerically for the Sabra shell model in Ref. [24]. We 
have carried out a similar exercise for the GOY shell model of fluid turbulence and 
found, unsurprisingly, that the universality of equal-time multiscaling exponents holds 
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Figure 3. (a) Representative plots of the normalised, third-order time-dependent 
structm'e function versus the dimensionless time t/tL from our numerical simulations 
of Model B for shell numbers n = 6 (uppermost) to n = 11 (lowermost). The inset 
shows a log- log plot of T|(A:„) versus A:„. The slope of this plot gives the dynamic 
sealing exponent = 1.398 ± 0.003, which agrees with Z2 from Fig. 2b. 
(b) A log- log plot of T^{kn) versus A:„. The slope of this plot gives the time-dependent 
sealing exponent Z4 — 1.402 ± 0.005. Our result from numerical simulations agree well 
with the analytical result Z2 ^ ~ 2 ~ S,, since wc have chosen ^ = 0.6. 



for the GOY model as well. Furthermore, we have also checked that a replacement of 
the viscous term z/q/c^ in equation (17) by the hyperviscous term i'aknikn/kd)°' does not 
affect these exponents [40]. 

We investigate now whether time- dependent velocity structure functions in decaying 
fluid turbulence show a factorisation similar to (32) and (37) for Models A and B. Let 
us normalise the order-p, time-dependent structure function by its value at t = 0: 

^ F^{kn,to,t)_ 

- Fp-{k^,to,Oy ^''^ 

to examine the dependence of Qp{n,t) on to- We give, in Figure (4a), representative 
plots of Q2{n = 9, t) and Q2{n = 5, t) versus t/t^ for 6 different time origins to; succesive 
values of to are separated from each other by OMl- To the extent that the symbols for 
different values of to overlap, our results show that Qp{n,t) is independent of to (so we 
have not included to as an argument of Qp), i.e., Fp{kn, to, t) factorises into a part that 
depends on to and another that does not [cf., (32) and (37) for the simple, linear Models 
A and B]. Clearly the dynamic-multiscaling exponents extracted from Qp{n,t) cannot 
depend on to- 

Furthermore, these exponents are the same as their counterparts for statistically 
steady, forced turbulence. We illustrate this in Fig. (4b) by comparing our numerical 
results for the normahsed, time-dependent structure function F{{kn,t)/S'^{kn), for 
statistically steady, forced (superscript /) turbulence, and Q2{n,t). [To obtain a 
statistical steady state for the GOY shell model we use the external force = 
{1 + i) X 5 X 10~^6n,i-] We have checked explicitly, for 1 < p < 6, that Qp{n,t) 
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and F^{kn,t)/ Sp{kn) agree within our error bars. Hence we propose the following 
factorisation that relates the time-dependent structure functions for decaying and 
statistically steady (superscript /) turbulence, 

F^{kn, to, t) = Ap{kn, to)F^{kn, t), (67) 

with all the to dependence on the right-hand side in the coefficient function Ap. If we 
now use the multifractal form for suggested in [6], we obtain the shell- model analogue 
of equation (64). 




Figure 4. (a) Representative plots of Qp{n,t) versus the dimensionless time t/t^, for 
p=2, n~9 (lower curve) and p = 4, n=5 (upper curve), and 6 different time-origins 
to in decaying fluid turbulence for Model C. Successive time origins are separated by 
0.5tL- The different symbols for the different sets of data for the different time-origins 
are indistinguishable; this is compelling numerical evidence for our proposed factorised 
form of the time-dependent structure functions. 

(b) Plots of QX{kn,t) and Fl{kn,t) as a function of t/tL to numerically test our 
proposed Ansatz (67) decaying turbulence in Model C. Results are shown only for 
shells n = 4 (uppermost curve), 6, 8, and 12 (lowermost curve) for clarity. 



We calculate the integral- and the derivative-time scales and their associated 



exponents 2;^'" and Zp^^. The counterparts of equations (49), for M = 1, and (50), 
for M = 2, for the GOY model are : 



Q;in,t)dt; 



dt^ 



n-1/2 



(68) 



(69) 



here t^ is the time at which Qp{n,t) = fi, with < yU < 1. In principle we should 
use /i = 0, i.e., = 00, but this is not possible in any numerical calculation since 
Qp cannot be obtained accurately for large t. We use /i = 0.6; and we have checked 
in representative cases that our results do not change for 0.3 < /i < 0.7. To compute 



D,ut 
P,2 



we use a centred, finite-difference, sixth-order scheme. Slopes of log-log plots 
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of Tp^^{n) and T^,^{n) versus /c„, (4 < n < 14) give us z^^ (Fig 5a) and 2 (Fig 5b), 
respectively. 

Our results for equal-time and dynamic-multiscaling exponents for the GOY 
model (for both types of initial conditions) are given in Tables 1 and 2, respectively. 
We compute the multiscahng exponents for equal-time and time-dependent structure 
functions for 50 different cases. Tables 1 and 2 list the means of these values and their 
standard deviations yield the error bars. By comparing columns 2 of Tables 1 and 2 
with column 2 of Table II in Ref. [6] we confirm, for the GOY model, the weak version 
of universality [24], i.e., the equal-time exponents are the same for both decaying and 
statistically steady turbulence. Furthermore, our exponents for the GOY shell model 
agree with those presented in Ref. [24] for the Sabra shell model. 
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Figure 5. (a) Representative plots of Qg(n,t) versus t/ti, for Model C; for clarity 
we show shell numbers n = 4 (uppermost), 6, 8 and 12 (lowermost). The inset shows 
rg'"(n) versus fc„ on a log- log scale. A linear fit yields the dynamic-multiscaling 
exponent Zg'" = 0.76 ± 0.02. 

(b) Representative plots of Q^{n,t) versus t/t^ for the same shell numbers as in 
(a). The inset shows 74 2" (n) versus kn on a log- log scale. A linear fit yields 
zfj" = 0.76 ±0.01. 



A comparison of the remaining columns of Tables 1 and 2 with their counterparts in 
Table II of Ref. [6] shows that this weak universality also applies to dynamic multiscaling 
exponents. Moreover, our direct numerical results for z^^ and z^^ (columns 4 and 6 
in Tables 1 and 2) agree with the bridge-relation values of these exponents (columns 3 
and 5 in Tables 1 and 2) that follow from equations (53) and (54) and (columns 2 of 
Tables 1 and 2). Note that the agreement between corresponding entries in Tables 1 and 
2 shows that our results are insensitive to the type of initial conditions we use. Finally, if 
we compare these Tables with Table II in Ref. [6], we find that our dynamic-multiscaling 
exponents for decaying, homogeneous, isotropic turbulence agree with their counterparts 
for the statistically steady case. Pictorial comparisons of the data in Tables 1 and 2 and 
the results of Ref. [6] are shown in Figs. (6a) and (6b) for integral-time (columns 
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4 in Tables 1 and 2) and derivative-time (columns 6 in Tables 1 and 2) exponents. 
Similarly Figs. (7a) and (7b) compare the dynamic-multiscaling exponents from our 
direct numerical simulations with the values predicted for them by the bridge relations 
(53) and (54) and the equal-time exponents O given in Tables 1 and 2. 




Figure 6. (a) Plots of z^'" , with error bars, versus p for Model C with data from 
[6], for statistically steady turbulence, and columns 4 of Tables 1 and 2, for decaying 
turbulence; these plots illustrate the agreement between the three sets of exponents, 
(b) We compare the derivative-time exponents from [6] and columns 6 of Tables 
1 and 2. As in (a), we find that the three sets of exponents agree, within error bars, 
with each other. 



0.8 




7 



Figure 7. (a) Plots of z^'" , with error bars, versus p for Model C with values obtained 
via the bridge relations (columns 3, Tables 1 and 2) and those obtained from our 
numerical simulations (columns 4, Tables 1 and 2). 

(b)Plots of -2^2"' with error bars, versus p for Model C with values obtained via the 
bridge relations (columns 5, Tables 1 and 2) and those obtained from our numerical 
simulations (columns 6, Tables 1 and 2). 
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orutii \p j 






Lu 




D,u 


1 


0.379 ± 0.008 


0.621 ± 0.008 


0.61 ± 0.03 


0.68 ± 0.01 


0.699 ± 0.008 


2 


0.711 ± 0.002 


0.66 ± 0.01 


0.68 ± 0.01 


0.716 ± 0.008 


0.723 ± 0.006 


3 


1.007 ± 0.003 


0.704 ± 0.005 


0.711 ± 0.001 


0.74 ± 0.01 


0.752 ± 0.005 


4 


1.279 ± 0.006 


0.728 ± 0.009 


0.734 ± 0.002 


0.76 ± 0.02 


0.76 ± 0.01 


5 


1.525 ± 0.009 


0.75 ± 0.02 


0.755 ± 0.002 


0.77 ± 0.02 


0.77 ± 0.02 


6 


1.74 ± 0.01 


0.78 ± 0.02 


0.78 ± 0.03 


0.77 ± 0.03 


0.78 ± 0.02 



Table 1. Our simulation results for Model C with Type I initial conditions. Order-p 
(column 1); equal-time exponents Cp (column 2); integral-scale dynamic-multiscaling 
exponent z^'^ (column 3) from the bridge relation (53) and the values of Cp in column 
2; Zp'^ from our calculation using time-dependent structure functions (column 4); the 
derivative-time exponents 2^2" (column 6) from the bridge relation (54) and the values 
of Cp in column 2; z^^ from our calculation using time-dependent structure function 
(column 7). The error estimates are obtained as described in the text. 



order (p) 


/■u 
Sp 


4'r[Eq.(53)] 


I,u 

^p,l 




D,u 
^p,2 


1 


0.380 ± 0.004 


0.620 ± 0.004 


0.60 ± 0.02 


0.690 ± 0.009 


0.687 ± 0.003 


2 


0.709 ± 0.003 


0.671 ± 0.007 


0.67 ± 0.03 


0.72 ± 0.01 


0.719 ± 0.005 


3 


1.000 ± 0.005 


0.709 ± 0.008 


0.707 ± 0.006 


0.74 ± 0.02 


0.743 ± 0.007 


4 


1.266 ± 0.008 


0.73 ± 0.01 


0.736 ± 0.008 


0.76 ± 0.03 


0.75 ± 0.01 


5 


1.51 ± 0.01 


0.75 ± 0.02 


0.752 ± 0.009 


0.77 ± 0.03 


0.77 ± 0.02 


6 


1.74 ± 0.02 


0.77 ± 0.03 


0.76 ± 0.02 


0.77 ± 0.03 


0.76 ± 0.02 



Table 2. Our simulation results for Model C with Type II initial conditions. Order-p 
(column 1); equal-time exponents Q (column 2); integral-scale dynamic-multiscaling 
exponent z^'" (column 3) from the bridge relation (53) and the values of Cp in column 
2; Zp" from our calculation using time-dependent structure functions (column 4); the 
derivative-time exponents 2:^2" (column 6) from the bridge relation (54) and the values 
of Cp in column 2; z^^ from our calculation using time-dependent structure function 
(column 7). The error estimates are obtained as described in the text. 

5.3. Model D 

From a numerical study of a passive-scalar shell model, with advecting velocities from 
the Sabra shell model, it was shown [37, 38] that the equal-time scaling exponents (p are 
universal: they are the same for decaying and statistically steady turbulence; and, in the 
latter case, they do not depend on the type of forcing. We find, not surprisingly, that 
this universality holds even when the advecting velocity field is a solution of the GOY 
model, i.e., for Model D: Table 3 column 2 shows the values of for 1 < p < Q, which 
we have obtained for decaying turbulence; these agree with the results of Refs. [7, 32] 
for statistically steady turbulence in Model D. Our equal-time exponents are also within 
error bars of their counterparts for the passive-scalar shell model of Refs. [37, 38]. We 
obtain from log-log plots such as Fig. (8a) for the modified, equal-time structure 



The Universality of Dynamic Multiscaling in Homogeneous, Isotropic Turbulence 25 

function (22) versus /c„; the slope of the hnear region 4 < n < 12 yields that is 
plotted versus p in Fig. (8b). 




Figure 8. (a) Representative plots of Sp(n) versus fc„ on a logarithmic scale for p = 
1 (uppermost curve) to 8 (lowermost curve) for Model D. 

(b) Plot of (^p, obtained from the linear, inertial region in (a), versus p. Our data 
points, shown as squares, are connected by a line; the error bars are smaller than the 
size of the symbol. The straight line corresponds to the Kolmogorov prediction of p/3. 



To analyse the time-dependent, passive-scalar structure functions we follow 
Section 5.2 for Model C: We obtain integral- and derivative-time scales from equations 
(49) and (50), for M = 1 and M = 2, respectively. In the integral in (49) we set the 
upper limit to t^, the time at which the normalised, time-dependent structure function 

is equal to /i, with < /i < 1. We use = 0.6, but we have checked in representative 
cases that our results remain unchanged if we use the range of values 0.3 < fi < 0.8. 
Slopes of log-log plots of Tp'^{n) versus kn yield z^'^. To extract the derivative-time scale, 
we use a centred, finite-difference, sixth-order scheme to obtain T^f{n) from which we 
get Zp 2 • Integral- and derivative-time multiscaling exponents are extracted from linear 
fits in the inertial range 4 < n < 10 as shown in the insets of the representative plots 
in Figs. (9a) and (9b) for, respectively, Ql{n,t) and Ql{n,t) versus t/t^. Finally, from 
the bridge relations (61) and our GOY-model, equal-time exponents C"i = —0.44 ±0.04 
and C2" = 0.709 ± 0.003 we obtain z^j = 0.56 ± 0.04 and z^f = 0.645 ± 0.003 in 
agreement with the values from our simulations listed in columns 3 and 4 of Table 
3. By comparing these columns with their counterparts in Table II of Ref. [7], we 
find agreement, within our error bars, between the dynamic- multiscaling exponents for 
Model D for both statistically steady and decaying turbulence. 
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Figure 9. (a) Representative plots of Q^r^{n^t) versus t/ii, for sliell numbers 6 
(uppermost) to 11 (lowermost) for Model D. The inset shows a log-log plot of Tg 2 (n) 
versus A linear fit yields = 0.562 ± 0.006. 

(b) Representative plots of (53(71, t) versus t/tL for the same shell numbers as in 
(a) . The inset shows T3 2 {n) versus fc„ on a log- log scale. A linear fit yields 
z^:^ = 0.646 ± 0.003. 



6. Conclusions 

We have systematised the study of the dynamic multiscahng of time- dependent structure 
functions in four Models (A-D) for passive-scalar (A,B, and D) and fluid (C) turbulence. 
By a suitable normalisation of these structure functions, we eliminate their dependence 
on the origin of time to at which we initiate our measurements. We have shown 
analytically that for the Kraichnan model of passive-scalar turbulence (Models A and B) 
the two-point time-dependent structure function can be factorised into two parts, one 
depending on the origin of time to and another that does not. This suggests a suitable 
normalisation of the structure functions by which we can eliminate their dependence 
on to- Surprisingly the same normalisation works for other models of fluid turbulence 
(Model C) and passive-scalar turbulence (Model D) as we have shown by extensive 
numerical simulations. Once this dependence on to has been factored out, the methods 
developed earlier [6, 7] for statistically steady turbulence yield linear bridge relations 
that connect dynamic-multiscaling and equal-time exponents. We show analytically, for 
Models A and B, and numerically, for Models B, C, and D, that these exponents and 
bridge relations are the same for statistically steady and decaying turbulence. Thus 
we have generalised the universality of equal-time exponents [24] and provided strong 
evidence for dynamic universality, i.e., dynamic-multiscaling exponents do not depend 
on whether the turbulence decays or is statistically steady for each of the Models A- 
D. We have shown elsewhere that these exponents do not depend on the dissipation 
mechanism by studying the GOY model with a hyperviscous term [40] . 

It is useful to distinguish our results from other studies of temporal dependences of 
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quantities in turbulence. Two such results are described below. 

The first example is the result Etot{t) ~ which holds for the total energy Etot at 
large times t in decaying turbulence [1]. This result holds for times that are longer than 
the time at which the integral length scale becomes comparable to the size of the system. 
[In decaying turbulence, the integral scale Lintit) = {J dkE{k,t)/k)/ {J dkE{k,t)) 
increases with time since the large-/c part of the energy spectrum E{k,t), at time t, 
gets depleted as t increases.] In all the results we report here, our shell-model analogue 
of Lint is well below the linear size of the system, i.e.. Lint ^ kQ^. Thus our results are 
not modified significantly by finite-size corrections, nor does the trivial decay of Etot, 
mentioned above, set in and mask the djTiamic multiscaling we have elucidated. 

The second example is the intermittency of velocity time increments studied 
in Ref. [41]. These studies calculate structure functions along a single Lagrangian 
trajectory. Spatial separations, of the sort used in defining equal-time Eulerian structure 
functions, are replaced by temporal separations along a Lagrangian trajectory. This is 
distinct from the spatiotemporal structure functions we consider and which are required 
for a full elucidation of dynamic multiscaling here (and analogous dynamic scaling in 
critical phenomena). 

We have described the way in which we have obtained error bars for the equal- 
time and dynamic-multiscaling exponents. These error bars only account for statistical 
errors but not the systematic errors associated with the values of n over which we fit 
inertial-range exponents. One can try to estimate such systematic errors by obtaining 
local slopes of log- log plots that yield these exponents. However, local slopes can be 
deceptive in shell models since the values of kn are separated by factors of 2. Instead, 
we can try to estimate these systematic errors by comparing the exponents we obtain 
by fitting over the ranges 3 < n < 15, 4 < < 14, and 5 < < 13. We have carried out 
such checks in representative cases for the dynamic-multiscaling exponents we report. 
The error bars we then obtain are about a factor of 4 larger than those shown in Tables 
1-3. 

We hope our work will stimulate experimental studies of dynamic multiscaling 
in turbulence. Recent advances in the experimental techniques of particle tracking 
in turbulent flows [42, 43] have made it possible to obtain accurate measurements of 
Lagrangian properties. To obtain the types of structure functions that we have described 
it will be necessary at least to track two Lagrangian trajectories of particles that are 
separated initially by a distance r. At the level of second-order structure functions this 
has been attempted in the direct numerical simulations of Ref. [19]. 
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Cp 


I.e 


D,e 

^p,2 


1 


0.342 ± 0.002 


0.522 ± 0.002 


0.632 ± 0.003 


2 


0.634 ± 0.003 


0.531 ± 0.004 


0.647 ± 0.003 


3 


0.873 ± 0.003 


0.553 ± 0.006 


0.646 ± 0.003 


4 


1.072 ± 0.004 


0.563 ± 0.003 


0.642 ± 0.005 


5 


1.245 ± 0.004 


0.562 ± 0.006 


0.643 ± 0.006 


6 


1.370 ± 0.006 


0.576 ± 0.006 


0.640 ± 0.005 



Table 3. Our simulation results for Model D. Order-p (column 1); equal-time 
exponents C,^ (column 2); integral-scale dynamic- multiscaling exponent Zj^\ (column 
3); the derivative-time exponents (column 4). The error estimates are obtained 
as described in the text. 
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